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Discrete Weierstrass Fourier Transform 
and Experiments 

Sheng Zhang, Brendan Harding 


Abstract 

We established a new method called Discrete Weierstrass Fourier Trans¬ 
form, a faster and more generalized Discrete Fourier Transform, to approx¬ 
imate discrete data. The theory of this method as well as some experi¬ 
ments are analyzed in this paper. In some examples, this method has a 
faster convergent speed than Discrete Fourier Transform. 


1 Introduction 

Barnsley, Harding, Vince, and Viswanathan introduced a notion of Weierstrass 
Fourier Series to approximate rough functions (see HD- 

Key idea of constructing Weierstrass Fourier Series is: 

1. Construct a linear operator on ZT(M); 

2. Using this linear operator, transform the classical Fourier basis to a new 
basis; 

3. By Gram-Schmidt process, get an orthonormal basis; 

4. Using this orthonormal basis, do what we can do in classical Fourier Analysis. 

In this paper, we are going to extend their theory by the following steps: 

1. Following the key idea in their paper, establish a new expression of Weier¬ 
strass Fourier Series, which is more suitable for subsequent discussions in this 
paper. 

2. Deduce the notion of Discrete Weierstrass Fourier Transform. 

3. Provide some numerical examples to test this transform. 


2 Weierstrass Fourier Series 


In classical Fourier Analysis, {e 27Tlkx }kez is a complete orthonormal basis for 
L 2 ([0,1]), and we can approximate functions with Fourier Series. Barnsley, 
Harding, Vince, and Viswanathan constructed another orthonormal basis for 
L 2 ([0,1]), and established a theory for approximating functions with Weierstrass 
Fourier Series. A key result in their paper is listed below (see Theorem 2.1). 
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Theorem 2.1 (Barnsley et al.) (see pQ, Corollary 3.1). Assume thatp G [1, oc] ; 
a, b G M, b ^ 0, and \a\ ^ \b\p. For any g G L P (M) ; there is a unique solution 
f G L P (R) to the equation 


f{x) - af(bx ) = s(z), (1) 

and the solution is given by the following series, which are absolutely convergent 
in L p (M); 

[ Em=o a m g(b m x) if\a\<\b\p 
f(x) = l 

{ -E~=i G) m 5(#r) i/N>|6|'. 

The following corollary is a direct conclusion from this theorem. 

Corollary 2.2. Assume that a G [0,1) and b = 2. If g G L°°(1R) ; then 0 /ms 
a unique solution f G L°°(M), and the solution is given by the following series, 
which are absolutely convergent in 

oo 

/(*)= ]T a m fl (2"V). (2) 

m =0 


Furthermore, f G L 2 ([0,1]). 


Proof. Let p = oc in Theorem |2.1[ and the first part of the corollary is proved. 
As / G L°°(R), / G L°°([0,1]). So, / G L 2 ([0,1]). □ 

In paper pQ, Barnsley et al. introduced Weierstrass Fourier Series by sub¬ 
stituting g(x) = sin kx and cos kx. However, another form shown below is more 
suitable for subsequent discussions in this paper. 

In 0. assume that a G [0, 1) a nd b = 2. For each k G Z, let g(x) = 


ek(x) = e 27rzkx . By Corollary 2.2, there exists a unique solution // ~(x) = 


« m e/c(2 m x) G L 2 ([0,1]). By normalizing {/&}, a normalized basis {e^} 
for T 2 ([0,1]) is established: 


if k = 0 


e/c = 


E^o a rn e 27rik - 2rnx if k ^ 0. 


By Gram-Schmidt process, an orthonormal basis {e^} for L 2 ([0,1]) is ob¬ 
tained: 


Ck = { ep 


= VT^ E^o a rn e 27rik - 2mx 


e k -ae k/2 _ _ 2\ m : 

- V 1 a f 2^m=0 a e 


m /D 2irik-2 rn x 


if k = 0 
if k is odd 
— ae* lkx if k is even and k ^ 0. 


Notice that when a = 0, for all k G Z. To sum up: 
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Theorem 2.3. Given a G [0,1), the set of functions {ek}kez is on orthonormal 
basis for L 2 ([ 0,1]). 

Weierstrass Fourier coefficients for a function / can be calculated as follows: 




if k = 0 


<4 = (/> 4 ) = < 


(/, 4 ) 


(f,e k )-a(f,e k/2 ) 

y/l—a 2 


= yr^E“=o« m ^-2- 

= (1 - a 2 )Em=0 amQ; fe-2™ 


if k is odd 

actk /2 if k is even and k ^ 0, 


where ak is the k th Fourier coefficient. Notice that when a = 0, ak = ak for all 
k G Z. Now, we can approximate any function in L 2 ([0,1]) using Weierstrass 
Fourier Series. 


Theorem 2.4. With the notations above, given a G [0,1), if f G L 2 ([0,1]) has 
Fourier expansion 

oo 

/(a;) = Y a k e k(x), 

k= — oo 

then it has Weierstrass Fourier expansion 


oo 

f(x) = Y & kek(x). 

k— — oo 


In particular, when a = 0, two expansions are the same. 

Proposition 2.5. With the notations above, given a G [0,1), if f G L 2 ([0,1]) 
is a real function, then the following statements hold: 

(1) ak = Oi-k f° r all k G Z; 

(2) Ck — e_k for all k G Z; 

(3) Y!k=-n < 444 ) is real and Ylk=-n < 444 ) = + 2 Re(YJk=i < 444 )) for 

all n G Z+. 


3 Discrete Weierstrass Fourier Transform 

With the notations above, given a G [0,1), for any fixed n G N+, let A be an 
n x n matrix such that 


Aij = 4(^),0 < hi <n-l. (3) 

Given n data points bo, 4, 4, * * * > 4-i G C, let b = ( bo , 4,4, • • • , 4-i) T - 

Discrete Weierstrass Fourier Transform(DWFT) is defined as the linear 
operator on C n : dwft(b) = A -1 b. Inverse Discrete Weierstrass Fourier 
Transform(IDWFT) is defined as the linear operator on C n : idwft(b ) = Ab. 
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Theorem 3.1. Given a G [0,1), dwft and idwft definded above are inverse 
linear operators on C n . In particular, when a = 0, dwft is classical Discrete 
Fourier Transform and idwft is classical Inverse Discrete Fourier Transform. 

Like DFT, DWFT is useful for data compression. Fix n G N+. Let k G N + 
and 1 < k < if n is odd, 1 < k < n ^ L if n is even. Assume 

b= {b 0 ,b 1 ,b 2 ,--- ,b n -i) T (4) 

is a set of data. Do DWFT on b and get 

c = dwft(b) = (c 0 ,ci,c 2 , • • • ,c„_i) T . (5) 


Change c k , c k + 1 , • • • , c n - k to 0 and get 

c'(k) = (c 0 ,ci, • • • ,c fe _i,0, • • • , 0, c n - k+ i, ■ ■ ■ ,c n -i) T . (6) 


Do IDWFT on c'(fc) and get 

b\k) = idwft(c'(k)) = (b' 0 , b[, b’ 2 , ■ ■ ■ , b' n _ 1 ) T . (7) 


When b = (bo, bi, & 2 , * * • , b n -i) T is real, b'(k) we got above is not necessarily 
real. However, we always want to get real approximation for real data. So, if 
the data are real, a last procedure that discarding the imaginary part of b'(k) 
should be taken, say 

b'\k)=Re(b\k)). (8) 


Say b"(k) is the approximation of b using k terms using DWFT. Nevertheless, 
the imaginary part of b'(k) is usually very small compared to its real part. 

Now, if b"(k) is very closed to b even if k is very small compared to n, we 
can store c'(fc), a, and n instead of data b. Then we can calculate matrix A as 
well as b"(k) when we need and use b"(k) as an approximation of b. In this way, 
plenty of space is saved. 

If DWFT and IDWFT are replaced by DFT and IDFT, then the approxi¬ 
mation using DFT is obtained. 

By Theorem 3.1, DWFT is a generalization of DFT. Since we can choose 
a suitable a for each set of data, DWFT will never be worse than DFT in 
approximation. 


4 Experiments 


In this section, some data sets are approximated using DFT and DWFT. With 
the notations in ©, <§, ©, 0 and given data vector b and a G [0,1), 
error vector for k terms is defined as 


E k = b"(k) - 6, 

and error function for / is defined as the function 

Kk) = m\\2, 


(9) 


( 10 ) 


which maps k to the 2-norm of E k , for k G {l,2,...,|-hl}. In particular, 
/i(§ + 1) = 0 by Theorem 


3.1 
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4.1 Data from real functions 

In this part, the number of data n = 1024 for all examples, and each data vector 
b = (/(0),/(^),/(|), • • • for some real-valued function / on [0,1]. 

4.1.1 Linear function 


f(x) = x — 0.5 and a = 0.5 



Figure 1.1.1 The figure of the error function for f{pc) = x — 0.5 with a = 0.5. The 
blue curve represents approximation using DFT, and the green curve represents 
approximation using DWFT. DWFT is always better. 
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Figure 1.1.2 Approximation of f{pc) = x — 0.5 using 10 terms. The red curve 
represents the original data, the blue curve represents approximation values 
using DFT, and the green curve represents approximation values using DWFT. 
See also Figure 1.1.3 and Figure 1.1.4. 



Figure 1.1.3 Approximation of f(x) — x — 0.5 using 50 terms. 
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Figure 1.1.4 Approximation of f(x) = x — 0.5 using 300 terms. 

4.1.2 Triangular function with low frequency and high frequency 

f{pc) = sinx + 0.01 sin 105x and a = 0.5 



Figure 1.2.1 The figure of the error function for f(x) = sin x +0.01 sin 105x with 
a = 0.5. The blue curve represents approximation using DFT, and the green 
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curve represents approximation using DWFT. DWFT is better at the first 510 
terms. 



Figure 1.2.2 Approximation of f(x) = sinx + 0.01 sin 105x using 10 terms. The 
red curve represents the original data, the blue curve represents approximation 
values using DFT, and the green curve represents approximation values using 
DWFT. See also Figure 1.2.3 and Figure 1.2.4. 
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Figure 1.2.3 Approximation of f(x) = sinx + 0.01 sin 105x using 50 terms. 



Figure 1.2.4 Approximation of /(x) = sinx + 0.01 sin 105x using 300 terms. 

4.1.3 Discontinuous function 

= * 6 ( 1 , 1 ] anda = 0 ' 5 
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Figure 1.3.1 The figure of the error function for f(x) = 


o x e [o, \] , 

1 wrth 

a = 0.5. The blue curve represents approximation using DFT, and the green 
curve represents approximation using DWFT. DWFT is not better. 



Figure 1.3.2 Approximation of f{x) — 


{? 


using 10 terms. The 


x e [0, |] 

* e (|, l] 

red curve represents the original data, the blue curve represents approximation 
values using DFT, and the green curve represents approximation values using 
DWFT. See also Figure 1.3.3 and Figure 1.3.4. 
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Figure 1.3.3 Approximation of f(x) 


0 x G [0, ~] 
1 x G (|, 1] 


using 50 terms. 
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Figure 1.3.4 Approximation of f(x) = 


0 x G [0, \\ 
l x e (|, l] 


using 300 terms. 
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4.1.4 Rough function 

f(x ) = Xlfclo 0-42^ cos(7r • 2 k x ) and a = 0.42 



Figure 1.4.1 The figure of the error function for f{x) = 0-42^ cos(7r • 2 k x) 

with a = 0.42. The blue curve represents approximation using DFT, and the 
green curve represents approximation using DWFT. DWFT is better at the first 
454 terms. 
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Figure 1.4.2 Approximation of f{pc) = XlfcLo 0-42^ eos( 7r • 2 k x ) using 10 terms. 
The red curve represents the original data, the blue curve represents approxi¬ 
mation values using DFT, and the green curve represents approximation values 
using DWFT. See also Figure 1.4.3 and Figure 1.4.4. 



Figure 1.4.3 Approximation of f{pc) = 0-42^ cos(7r • 2 k x) using 50 terms. 



Figure 1.4.4 Approximation of f{x) = YkL o 0.42 fe cos(7r • 2 k x) using 300 terms. 
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4.2 Discrete data 


In this part, data are practical data from some websites. Given data vector b 
and a E [0,1), error vector and error function for the data are defined in § 
and (fTol) . 


4.2.1 Stock price of Commonwealth Bank of Australia 

Daily open prices of Commonwealth Bank of Australia from Sep 30, 2010. To¬ 
tally 1024 data points. Source from https://au.finance.yahoo.com/. 



Figure 2.1.1 The figure of the error function for 1024 daily open prices of Com¬ 
monwealth Bank of Australia with a = 0.3. The blue curve represents approxi¬ 
mation using DFT, and the green curve represents approximation using DWFT. 
DWFT is better at the first 88 terms. 
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Figure 2.1.2 Approximation of 1024 daily open prices of Commonwealth Bank 
of Australia using 10 terms. The red curve represents the original data, the 
blue curve represents approximation values using DFT, and the green curve 
represents approximation values using DWFT. See also Figure 2.1.3 and Figure 
2.1.4. 



Figure 2.1.3 Approximation of 1024 daily open prices of Commonwealth Bank 
of Australia using 50 terms. 
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Figure 2.1.4 Approximation of 1024 daily open prices of Commonwealth Bank 
of Australia using 300 terms. 

4.2.2 Water level of Alameda in California 

Water level of Alameda in California per hour from May 29, 2014. Totally 1024 
data points. Source from http://tidesandcurrents.noaa.gov/. 



16 












Figure 2.2.1 The figure of the error function for 1024 water levels of Alameda in 
California per hour from May 29, 2014 with a = 0.3. The blue curve represents 
approximation using DFT, and the green curve represents approximation using 
DWFT. DWFT is not better. 
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Figure 2.2.2 Approximation of 1024 water levels of Alameda in California using 
10 terms. The red curve represents the original data, the blue curve represents 
approximation values using DFT, and the green curve represents approximation 
values using DWFT. See also Figure 2.2.3 and Figure 2.2.4. 
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Figure 2.2.3 Approximation of 1024 water levels of Alameda in California using 
50 terms. 
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Figure 2.2.4 Approximation of 1024 water levels of Alameda in California using 
300 terms. 
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4.3 Comments 

4.3.1 


Figure 1.1.2 and Figure 1.1.3 in Section 4.1.1 coincide with Figure 1, Figure 2, 
Figure 3 and Figure 4 of [I], which means that the discrete algorithm in this 
paper gives the same result as the continuous algorithm when the number of 
the data is much bigger than the terms of approximation. 

4.3.2 

In the figures above we can see that at some of the data points, DWFT approx¬ 
imates better, while at some other data points DFT is better. For example, in 
Figure 1.1.3, at the data points near 0.0 or 1.0, DWFT is better while at the 
data points near 0.5, it is significantly worse, but overall, it approximates the 
data better in the norm defined above. 


4.3.3 


We can also see a characteristic that approximation values using DWFT shakes 
heavily at the points near |, less heavily at the points near |, |, and less heavily 
at the points near ^|, f, -, 


4.3.4 

In the above examples, DWFT is better than DFT for some smooth continuous 
functions. For some other smooth functions such as /(x) = x 2 with a = 0.5 and 
/(x) = e x with a = 0.5, DWFT is also better than DFT. But for /(x) = x(x —1) 
with a = 0.5, DWFT behaves very bad. 

4.3.5 

For triangle functions with low frequency and high frequency such as the exam¬ 
ple in Section 4.1.2, /(x) = sinx + 0.01 cos 100x + 0.01 cos200x with a = 0.5 and 
/(x) = 10 sin 0.lx — cos x + 0.01 sin lOOx + 0.02 cos 200x with a = 0.5, DWFT is 
better than DFT. 

However, make sure that the low frequencies are low enough with coefficients 
big enough compared to the high frequencies, and the low frequencies are not 
counteracted each other. For instance, DWFT doesn’t work better than DFT for 
the function /(x) = 10 sin 0.lx — sin x +0.01 sin 100x + 0.02 cos 200x with a = 0.5 
because the low frequency terms “10 sin 0.lx” and “—sinx” are counteracted 
each other. 


4.3.6 

For discontinuous functions, DWFT doesn’t seem better than DFT. As we can 
see in Section 4.1.3, approximation values using DWFT for discontinuous func¬ 
tion shake heavily. 
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4.3.7 


As we can see in Section 4.1.4, DWFT has some advantages in approximating 
self-similar rough functions. 


4.3.8 

As we can see in Section 4.2, DWFT seems not good for approximating practical 
data, because practical data are not perfectly self-similar. A small shake in the 
data can induce shakes everywhere in the approximation values using DWFT. 
In this case, the approximation values using DWFT look much rougher than 
the data. 


4.3.9 

If some data set can be approximated better using DWFT, we have a better 
way to compress the data. However, the calculation of DWFT is much more 
complicated than DFT. We might wish to find a fast way to calculate DWFT. 
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